Second-order nonlinear optical effects of spin currents 
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Abstract 

A pure spin current formed by opposite spins moving in opposite directions is a rank-2 axial tensor 
which breaks the inversion symmetry. Thus a spin current has a second-order optical susceptibility, with 
unique polarization-dependence determined by the symmetry properties of the current. In particular, a 
longitudinal spin current, in which the spin polarization directions are parallel or anti-parallel to the moving 
directions, being a chiral quantity, leads to a chiral sum-frequency effect. Microscopic calculations based 
on the eight-band model of a III-V compound semiconductor confirm the symmetry analysis and show that 
the susceptibility is quite measurable under realistic conditions. The second-order nonlinear optical effects 
may be used for in-situ and non-destructive detection of spin currents, as a standard spectroscopy tool in 
research of spintronics. 

PACS numbers: 72.25.Dc, 42.65.An, 78.20.Ls 
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spin currents, which carry information via spins in heu of charges, play a key role in spin- 
tronics (1112]]. Pure spin currents also signify the occurrence of some novel spin-related quantum 
phenomena such as the spin Hall effect OUTO]!. the quantum spin Hall effect and topological insula- 
tors ifTTUTSl . Spin currents were previously observed via spin accumulation at stopping edges 0- 
[9l [T6l or conversion to electrical signals ifTOl [T7l - [T9l . Direct and non-destructive measurement 
of pure spin currents where and while they flow [[20l [2T]] is highly desired, but is very difficult 
because a pure spin current bears neither net charge current nor net magnetization. Noticing that 
a longitudinal spin current in which the spins point parallel or anti-parallel to the current is a chi- 
ral quantity, we envisaged that it can be probed by the chiral sum-frequency optical spectroscopy 
which was recently developed to detect molecular chirality [l22l - [24ll . By symmetry analysis in gen- 
eral cases and microscopic calculations in realistic models, we discovered that a pure spin current 
has sizable second-order optical susceptibility. This finding lays the foundation of direct, non- 
destructive measurement of spin currents by standard optical spectroscopy, facilitating application 
of spintronics [[1113 and research on spin-related quantum phenomena Ol-fTSll. 

As a basic principle of nature, a physical object is measurable only when it breaks certain fun- 
damental symmetries. Indeed, the probe must break the same symmetries as the object does, since 
the whole coupled system of an object and a probe has the fundamental symmetries. For example, 
in an Ampere meter, a "pure" charge current, which breaks the time reversal symmetry, is cou- 
pled to a microscopic current inside a magnet. Such symmetry consideration led to a scheme of 
detecting a pure spin current by a "photon spin current" carried by a polarized light beam Il20ll . 
A recent experiment [l2T1l showing coupling between a spin current and a spin wave is a remark- 
able demonstration of the symmetry principle of measurement. Though as a direct probe of spin 
currents, the spin-wave technique [l2T1l still requires special design and fabrication of magnetic 
nanostructures and the "photon spin current" probe [20] is limited by weak interaction since it 
involves the tiny light momentum, these previous works paved the way of searching methods of 
direct and non-destructive measurement of pure spin currents using symmetry analysis. 

Spin currents have peculiar symmetry properties owing to the characteristics of spins. Unlike 
a charge which is a scalar, a spin is a vector pointing to a certain direction. Physically, a spin is 
like a tiny magnet resulting from a quantized amount of current circulating about the spin direc- 
tion [Fig.[T](a)]. Such physical nature makes a spin an unusual vector, namely, an axial vector. As 
illustrated in Fig.[T](a), a spin reverses inside a parallel mirror and is unchanged inside a perpendic- 
ular mirror, in opposite to a polar vector. Spin-ups and spin-downs moving in opposite directions 
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FIG. 1. (Color online) Symmetry analysis for sum-frequency effects of spin currents, (a) A spin under 
mirror reflections, (b) A longitudinal spin current under mirror reflections. The arrows indicate the moving 
directions of the spins, (c) Chiral sum-frequency processes in chiral systems, (d) A transverse spin current 
under mirror reflections. 

v^ith the same velocity make up a pure spin current w^ithout net spin polarization or magnetization. 
When the spins are parallel or anti-parallel to the moving direction, the spin current is a longitudi- 
nal one. A longitudinal spin current has a special symmetry property - chirality. An object, such 
as a hand or a helix, is chiral if it cannot be made identical to its mirror image by translation and 
rotation. The chirality of a longitudinal spin current is illustrated in Fig.[T](b): If a spin's micro- 
scopic current circulates its moving direction left-handedly, the mirror image does right-handedly, 
and vice versa. 

Noticing its chirality, v^e conceived the idea of measuring a longitudinal spin current using 
the chiral sum-frequency optical spectroscopy, v^hich was recently developed as a standard tool 
to study molecular chirality [l22l - [24ll . In chiral sum-frequency, two input optical fields Fi and 
F2 (with frequencies oji and (jl>2, respectively) and the induced polarization field P at frequency 
oji + (jj2 form a left- or right-hand system. Fig.[T](c) show^s how^ a chiral sum-frequency process 
occurs in a chiral system. Considering a right-hand helix, a charge at position A vv^ill be driven 
to point B by an electric field Fi w^hich is along the X-axis, and then to point C by F2 w^hich 
is along the F-axis. The confinement of the helix leads to a net displacement along the Z-axis. 
Thus the two input fields and the induced polarization (Fi, F2, P) form a right-hand system. If the 
order of input fields is reversed (F2 applies before Fi), the charge w^ould follow^ a trajectory like 
D ^ E ^ F, resulting in a polarization along the -Z-axis, and (F2, Fi, P) still form a right-hand 
system. Similarly, the sum-frequency in a left-hand helix is a left-hand chiral process. A sum- 
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frequency process is characterized by a second-order susceptibility x via P (oji + (02) = X ' 
Fi ((jJi)F2 ((JJ2)' In the above example of helix, the susceptibility may be written as a form of 
three dyadic vectors, = A (ZYX - ZXY), i.e., a rank-3 tensor. Thus the chiral sum-frequency 
susceptibility provides a measurement of the chirality of a physical object. If otherwise measured 
in linear optics, the effect of the molecular chirality relies on the small magnetic moment of the 
molecules, and in turn on the small photon momentum of the probe light, similar to the case of 
linear optical effects of spin currents GOll . 

For a systematic symmetry analysis, we consider a spin current with both longitudinal and 
transverse components. We define the Z-axis as the current direction and the X-axis as the spin 
direction of the transverse component. The spin current can be written as a rank-2 tensor J = 
/xXZ+/zZZ, in a form of dyadic vectors, in which the left/right vector is the spin/current direction 
and Jz/x is the longitudinal/transverse amplitude. Above all, the spin current breaks the inversion 
symmetry, satisfying the symmetry properties required by a second-order optical process [|25ll . 

In general, the sum-frequency susceptibility tensor has 27 independent ttrms^x^^^ = A'xxxXXX+ 
XxxY^^^ + • • • +A'zzzZZZ, but the symmetry properties of a spin current will set many terms to be 
zero or non-independent [|25l|. For a longitudinal spin current, only the chiral terms are non-zero. 
In a non-chiral term, at least one of the three directions X, Y and Z appears even times (twice or 
zero times). Consider ;i;xxFXXY for example. Under reflection by the Y-Z plane, the longitudinal 
spin current is reversed, hutxxxY^^^ is unchanged, so this term must be zero. Similar arguments 
apply to other non-chiral terms. Also, the susceptibility must be anti- symmetric under reflection 
by any plane parallel to the Z-axis. With these constraints, the sum-frequency susceptibility of a 
longitudinal spin current can be written as 

xf^ = Jz[ai (XYZ-YXZ) 

+a2 (YZX - XZY) + a3 (ZXY - ZYX) ], (1) 

with only three independent parameters. As for a transverse spin current /^XZ, it changes its sign 
under reflection by the X-Z plane but is invariant under reflection by the X-Y or Y-Z plane [see 
Fig. [T](d)], each non-zero term in the susceptibility must contain odd times of Y and even times of 
Z or X, so 

= Jx(xiXXY + JC2XYX + JC3YXX 

+Z1ZZY + Z2ZYZ + Z3YZZ + yYYY), (2) 
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FIG. 2. (Color online) Models for microscopic calculation of the sum-frequency susceptibility, (a) The 
full eight-band model and the electron spin distribution for a pure spin current in a semiconductor, (b) 
and (c) Selection rules and relative dipole moments from the spin-3/2 and spin- 1/2 valence bands to the 
conduction band, respectively, (d) A simplified model with the HH-LH splitting neglected. The spin states 
and selection rules for inter-band transitions are independent of the momentum. The transition energies to 
the Fermi surface from different valence bands are indicated. 

with seven independent parameters. The unique polarization-dependence of the second-order sus- 
ceptibility of a spin current can be used to distinguish its transverse and longitudinal components, 
and also to single out the spin-current signature from the effects of the material background or a 
charge current [ j26] |. 

To determine the independent parameters of the susceptibility in Eqns. ([T]) and (|2]), wt per- 
formed microscopic calculation for a pure spin current in a bulk GaAs, using the standard per- 
turbation theory [[25l[23 w^ith an eight-band model Il28ll . We assumed that the pure spin current 
result from a non-equilibrium distribution of electrons in the conduction band, with a small por- 
tion of electrons near the Fermi surface having opposite spin polarizations for opposite velocities 
[Fig. [2] (a)] , under conditions similar to those in Ref. [jTl. The optical interaction includes the 
inter-band transitions between the valence bands and the conduction band and the intra-band ac- 
celeration of electrons and holes. To avoid real absorption of light, the light frequencies were 
chosen such that the sum frequency is below the band gap. For the sake of simplicity, we ne- 
glected the anisotropy of the valence bands. We also adopted the free-particle approximation. 
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which is justified since the Coulomb interaction is largely screened in the n-doped material. These 
approximations, according to the symmetry analysis, would only quantitatively modify the results. 
The spin splitting due to the bulk inversion asymmetry of the material (the Dresselhaus eff'ect) is 
as small as 0.01 meV for the doping level considered (3 x 10^^ cm"^), and therefore was neglected 
in the calculation. The bulk inversion asymmetry would cause a background second-order sus- 
ceptibility, which is indeed strong but can be well separated from the spin-current eff'ect by AC 
modulation of the current and phase-locking detection. Two representative results of the calcu- 
lated susceptibility spectra are shown in Figs. [3] (a) and (b). The other terms of the susceptibility 
tensor (not shown) have similar frequency-dependence and comparable amplitudes. As a specific 
example, a transverse spin current 20 nA/yum"^ has a susceptibility ;t/7zz ~ 0.40 x 10"^ esu (or 
0.17 X 10"^^ m/V in SI units) for input frequencies oji = 100 meV and (jJ2 = 1,400 meV, or 
17. X 10"^^ esu for oji = (jj2 = 750 meV (corresponding to the second harmonics generation). 

To better understand the microscopic mechanism of the sum- frequency eff'ect of a spin current, 
we simplify the model by neglecting the splitting between the heavy hole (HH) band and the light 
hole (LH) band. Under this approximation, the HH and LH bands form a spin-3/2 band with 4- 
fold degeneracy. The split-off' (SO) band and the conduction band have spin- 1/2. In this simplified 
model, the spin states and the selection rules for inter-band transitions are separated from the 
momentum [Figs. [2] (b) and (c)]. 

Let us first consider a single electron with momentum k and spin polarization Sk [Fig.[2](d)]. We 
set up a coordinate system (ei, 62, 63) so that Sk = 63 (/+ - /_) /2 with /+/_ denoting the population 
at the spin-up/down state. The angular momentum conservation requires that a light with circular 
polarization ei ± ie2 couples only to the transitions |j, m) |l/2, m ± 1), where j = 3/2 or 1/2 
is the spin of a valence band and m = -j, j + 1, . . ., or j is the component along the e3-axis. The 
relative dipole moments of the relevant inter-band transitions are indicated in Figs.|2](b) and (c). 
To simplify the discussion, we set the input frequency (jl>2 to be near resonant with the band gap 
and much greater than cji, so that the inter-band transitions and the intra-band driving are mostly 
caused by F2 exp {-ico2h) and Fi exp {-icoiti), respectively. 

The probability amplitude of a certain inter-band transition is determined by the inner product 
of the dipole moment and the optical field. For example, the transition |3/2, -3/2) |l/2, 1/2) 
generated between t2 and t2 + dt2 has a probability amplitude dG2 = /(I - /-) [dl^l (ci + 1^2) * 
F2 exp(-/6(;2^2)^^2, where d^y is the inter-band dipole, and the factor (1 - /_) accounts for the 
Pauli blocking. After the excitation, the probability amplitude oscillates in time with frequency 
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£'3/2(k), leading to the optical polarization (ei + i^2)\dcvl ^/ij e~'^^^^^^^^^~^^^dG2 at time t, where 
£'3/2(k) = k^/(2me) + k^/(2m3/2) is the transition energy of a pair of electron an hole with mass 
rrie and m3/2, respectively. The radiation has the same circular polarization as the input because of 
the angular momentum conservation. Summation over all possible transitions and integration over 
time give the linear optical response to the field F2 as 




-iE3/2(k)(t-t2) 



X 2 (1 - A) (ei + ie2) (ei + /e2)* • F2^-^'"^^^ J^2. (3) 
+ 

Thus P^^^ oc (6162 - 6261) • F2 = F2 X Sk, which has a transparent physical meaning: The linear 
polarization of the output field is related to that of the input one by a rotation about the spin, 
essentially a Faraday rotation due to the spin acting as a magnet. When the eff'ect of the intra- 
band driving by Fi is included, the momentum k should be replaced with the accelerated one 
k-r = k - eFi txp(-i(jjiti)dti at time r, and the phase £'3/2(k)(^ - ^2) accumulated from the 
creation time t2 to the recombination time t should be replaced with £"3/2 ^k^] dr. By expansion 
to the linear order of Fi, we have ^ - 2ek • Fi j^^ txp(-i(jjiti)dti, so the second-order optical 
response can be written as P oc F2 x Sk^Vk • Fi, where Vk = k/m^ is the velocity of the electron 
with momentum k. The physical meaning of ^Vk • Fi is obviously the power done by the field to 
the electron. ^SkVk is just the spin current tensor contributed by the electron. 

For a distribution of electrons, the summation over the momentum space gives the sum- 
frequency response as P = ^F2 x (J • Fi), with 

^ ^ / g, + 2 \' (2/3) |rfevl'(l+m,/m3/2) 
\ 3 I {aJi-\-aj2- £3/2) {(^2 - E3/2) oji 

- 1^3/2, m3/2 £1/2, mi/2 j, (4) 

derived by Fourier transformation of Eq. ([3]) including the intra-band driving and the contribution 
of the SO band, where the factor containing the material dielectric constant Sr takes into account 
the diff'erence between the macroscopic external field and the microscopic local field Il29l , rnj 
denotes the mass of the spin- j hole band, and Ej is the transition energy from the spin- j band to 
the Fermi surface [see Fig. [2](d)]. The constants in Eqns. ([T]) and ([2]) are such that ai = -Z2 = 
Z3 = ^ and others= 0. With the HH-LH splitting neglected, the sum-frequency susceptibility has 
a compact form with only one independent parameter. This feature is due to the separation of 
the spin and motion degrees of freedom of the electrons and holes. When the HH-LH splitting 
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FIG. 3. (Color online) Representative results of the sum- frequency susceptibility. (^)xxyz due to a longi- 
tudinal spin current, and (b) Xyzz due to a transverse spin current, as functions of the optical frequencies. 
Parameters are chosen similar to those in Ref. Q: The band gap is 1519 meV, the HH-SO splitting is 
341 meV, the doping concentration is 3 x 10^^ cm"^, the effective mass (in units of free electron mass) of 
the HH, LH, SO, and conduction bands is in turn 0.45, 0.082, 0.15, and 0.067, the dipole dcv = 6.7 eA, the 
dielectric constant Sr = 10.6, and the spin current Jx = Jz = 20 nA/^m^. 

is non-zero, the spin quantization direction and therefore the optical selection rules depend on 
the momentum and vary w^ith acceleration of the particles. This leads to the general form of 
susceptibility in Eqns. ([T]) and ([2]), wiih the extra terms proportional to the HH-LH splitting. 

In summary, with systematic symmetry analysis in general cases and microscopic calculation 
under realistic conditions, wt have shov^n that a pure spin current has a sizable sum-frequency 
susceptibility. In particular, a longitudinal spin current has a chiral sum-frequency effect. The 
current results can be straightforv^ardly extended to other second-order optical spectroscopy such 
as difference-frequency and three- w^ave mixing ll25ll . Thus the standard nonlinear optical spec- 
troscopy makes up a toolbox for research of spintronics. With universality of the method guar- 
anteed by the symmetry principle and v^ithout requirements of resonance conditions or special 
structure design and fabrication, the nonlinear optical spectroscopy can be applied to study a v^ide 
range of spin-related quantum phenomena such as the quantum spin Hall effect and topological in- 
sulators ifTTUTSl . A vv^ealth of physics connecting spins and photons and technologies synthesizing 
spintronics and photonics are still to be explored. 
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